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ABSTRACT 

We present a formalism to study tidal interactions in white dwarf binaries in the limiting case of quasi-static 
tides, in which the tidal forcing frequencies are small compared to the inverse of the white dwarf's dynamical 
time scale. The formalism is valid for arbitrary orbital eccentricities and therefore applicable to white dwarf 
binaries in the Galactic disk as well as globular clusters. In the quasi-static limit, the total perturbation of the 
gravitational potential shows a phase shift with respect to the position of the companion, the magnitude of which 
is determined primarily by the efficiency of energy dissipation through convective damping. We determine rates 
of secular evolution of the orbital elements and white dwarf rotational angular velocity for a Q3Mq helium 
white dwarf in binaries with orbital frequencies in the Laser Interferometer Space Antenna (LISA) gravitational 
wave frequency band and companion masses ranging from O.3M0 to IQ^Mq. The resulting tidal evolution time 
scales for the orbital semi-major axis are longer than a Hubble time, so that convective damping of quasi-static 
tides need not be considered in the construction of gravitational wave templates of white dwarf binaries in 
the LISA band. Spin-up of the white dwarf, on the other hand, can occur on time scales of less than lOMyr, 
provided that the white dwarf is initially rotating with a frequency much smaller than the orbital frequency. 
For semi-detached white dwarf binaries spin-up can occur on time scales of less than 1 Myr Nevertheless, the 
time scales remain longer than the orbital inspiral time scales due to gravitational radiation, so that the degree 
of asynchronism in these binaries increases. As a consequence, tidal forcing eventually occurs at forcing 
frequencies beyond the quasi-static tide approximation. For the shortest period binaries, energy dissipation is 
therefore expected to take place through dynamic tides and resonantly excited ^-modes. 
Subject headings: Stars: Binaries: Close, Stars: White Dwarfs, Stars: Oscillations 



1. INTRODUCTION 

White dwarfs are the most common endpoint of stellar evo- 
lution in galaxies and dense stellar systems. During the past 
decade, ongoing optical surveys such as the Sloan Digital Sky 
Survey (SDSS) and the ESO Supernova la Progenitor surveY 
(SPY) have provided a plethora of new white dwarf binaries 
suitable to study poorly understood binary evolution phases 
such as common envelope evolution and type la supernovae. 
Binaries consisting of two white dwarfs are furthermore the 
single most abundant and only guaranteed sources of gravita- 
tional wave radiation in the Galaxy for the Laser Interferom- 
eter Space Antenna, LISA (Bender et al. 1998). 

Despite the high abundance of white dwarf binaries in 
present-day and planned astrophysical surveys, the study of 
their orbital evolution has remained limited to models in 
which the white dwarfs are treated as point masses, in so far 
as the calculation of the gravitational potential is concerned. 
Gravitational wave forms of double white dwarf binaries, for 
instance, currently only account for orbital frequency changes 
driven by gravitational radiation. As a significant fraction of 
these systems spirals in to periods as short as 5 minutes, tidal 
effects can, in principle, have a significant impact on the or- 
bital evolution and thus the gravitational wave frequency evo- 
lution. The impact depends strongly on the strength and the 
nature of the tidal energy dissipation mechanism, which, for 
white dwarfs, remains uncertain. Racine, Phinney, & Arras 
(2007) recently also proposed a non-dissipative tidal evolu- 
tion mechanism in which angular momentum exchange be- 
tween the star and the orbit is driven by resonant excitation of 
Rossby modes. However, as commented by the authors, non- 
linear effects likely limit the lifetime of the resonance and thus 
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the tidal evolution mechanism. 

Campbell (1984) studied dissipative tides in white dwarf bi- 
naries assuming the orbital frequency of the binary to be much 
smaller than the frequencies of the white dwarf's free modes 
of oscillation and limiting his study to circular binaries. The 
author considered energy dissipation through perturbations of 
the radiative energy flux, and used a perturbation technique to 
calculate the tidal velocity field in a non-rotating white dwarf. 
Campbell found that, in circular binaries, the synchronization 
time scale of a white dwarf can be comparable to the white 
dwarf's lifetime, provided that the initial degree of asynchro- 
nism is sufficiently large. However, in his estimate for the 
tidal synchronization time scale, he used a rather high white 
dwarf luminosity of 0.03 L©, which decreases the time scale 
by several orders of magnitude compared to the tidal synchro- 
nization time scale of a more conventional ~ Lq white 
dwarf (see his Eq. 46). Other authors studying the impact 
of tidal effects on the evolution of white dwarf binaries either 
assumed tidal dissipation to be strong enough to maintain syn- 
chronism at all times (Webbink & Iben 1987; Mochkovitch & 
Livio 1989; Iben, Tutukov, & Fedorova 1998) or parameter- 
ized tidal dissipation by means of an ad-hoc tidal synchroniza- 
tion time scale (Marsh, Nelemans, & Steeghs 2004; Gokhale, 
Peng, & Frank 2007). 

In contrast, tidal evolution theories for binaries with non- 
degenerate component stars are well developed and able to 
pass stringent observational tests such as measured orbital de- 
cay rates of high-mass X-ray binaries (Belczynski et al. 2008) 
and circularization cut-off periods in young open cluster bi- 
naries (Witte & Savonije 2002, Zahn 2008). However, for the 
latter test, some debate still exists on the agreement between 
theory and observations (Mathieu, Meibom, & Dolan 2004; 
vid^Mbfltsa'aaeiWisaoieu 2005). According to our current under- 
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standing, non-degenerate stars with convective envelopes dis- 
sipate tidal energy primarily through convective damping of 
quasi-static tides, while non-degenerate stars with radiative 
envelopes dissipate tidal energy primarily through radiative 
damping of dynamic tides (Zahn 1975, 1977). Both dissi- 
pation mechanisms may be significantly enhanced by reso- 
nances between tidally forced oscillations and free oscillation 
modes of the component stars, especially in binaries with ec- 
centric orbits (Witte & Savonije 1999, 2001, 2002). 

Our aim in this paper is to initiate a systematic investigation 
of tidal dissipation mechanisms operating in white dwarfs by 
building on the success of tidal evolution for non-degenerate 
stars. Here we explore the effectiveness of convective damp- 
ing of quasi-static tides as a tidal energy dissipation mecha- 
nism in white dwarf binaries. In § 2 and § 3, we outline the ba- 
sic assumptions and introduce the system of differential equa- 
tions governing tidally forced oscillations in close binaries. In 
§ 4, we use a perturbation method to derive approximate so- 
lutions to the system of differential equations appropriate for 
quasi-static tides. In § 5 and § 6, we present the equations gov- 
erning the secular evolution of the orbital elements and white 
dwarf rotation rates. In § 7, we calculate the quasi-static tidal 
distortion and orbital evolution time scales for a 0.3 He 
white dwarf model, and compare the orbital evolution time 
scales with those due to gravitational radiation. The final sec- 
tion is devoted to concluding remarks. 

2. THE TIDE-GENERATING POTENTIAL 

We consider a close binary system of stars revolving around 
one another in a Keplerian orbit with period Poib, semi-major 
axis fl, and eccentricity e. The first star, with mass M\, ra- 
dius R\, and luminosity Li, rotates uniformly around an axis 
perpendicular to the orbital plane with angular velocity ili in 
the sense of the orbital motion. The second star, hereafter re- 
ferred to as the companion, has mass and is considered to 
be a point mass. The rotational angular velocity 17i of star 1 is 
assumed to be small enough for the Coriolis force and the cen- 
trifugal force to be negligible, so that the tides raised by the 
companion can be treated as forced perturbations of a non- 
rotating, spherically symmetric, equilibrium star. 

The tidal force exerted by the companion is derived from 
the tide-generating potential £tide t) as 



with 



^lide = -Etide 



Etide =1 — ■ 

a j M\ 



(1) 



(2) 



The quantity £tide is a small dimensionless parameter corre- 
sponding to the ratio of the tidal force to the gravity at the 
star's equator. We express the tide-generating potential in 
terms of spherical coordinates r = (r,6,(p) with respect to an 
orthogonal frame of reference that is corotating with the star 
The polar angle 6 is measured from the rotational angular ve- 
locity vector, while the azimuthal angle (f) is measured in the 
orbital plane in the sense of the orbital motion. At f = 0, the 
= direction coincides with the direction from the star's 
mass center to the periastron of the binary orbit'. 

' For binaries with circular orbits, the </> = direction coincides with the 
direction from the star's mass center to the ascending node of the binary orbit 
at r = 0. 



As is customary, we expand the tide-generating potential in 
Fourier series as 

Etide W (r, t) = -Etide X! X! X! 



=2 m=-ek= 



— ) Yl"{9,4,)exp[i{a„at-knT) 



(3) 



(e.g. Polfliet & Smeyers 1990). In this expansion, G is 
the Newtonian constant of gravitation, Y"'{6,(j)) an unnormal- 
ized spherical harmonic of degree £ and azimuthal number m, 
o'm.k = kn + mili a forcing angular frequency with respect to 
the corotating frame of reference, « = 2 tt/ the mean mo- 
tion, and T a time of periastron passage. The factors Q,m,jfc are 
Fourier coefficients defined as 

ce.m.k = ,„ . I ,, Pf (0) — 



{£+\m\y. 



(1- 



£-1/2 



1 



X— / (1+ecosv)^ ' cos(^M+mv)^^v, 



(4) 



where v is the true anomaly, M = n(t- r) the mean anomaly, 
and P'"{x) an associated Legendre polynomial of the first 
kind. The coefficients ,„ji obey the property of symmetry 

C( -,„ -k = ci,,„± and, since f j,"''(0) = for odd values of (.+ \m\, 
are equal to zero for odd values of £ + |m|. From the binomial 
theorem it furthermore follows that .,„.o = for m = ±£. For a 
given orbital eccentricity and sufficiently large Fourier indices 
k, the non-zero coefficients .„, jt decrease with increasing val- 
ues of k, though the decrease is slower for higher orbital ec- 
centricities^ (e.g. Smeyers et al. 1998, Willems 2003). There 
is thus only a finite number of non-negligible terms contribut- 
ing to the expansion of the tide- generating potential, and the 
number of non-negligible terms increases with increasing or- 
bital eccentricities. 

The expansion of the tide-generating potential in Fourier 
series introduces an infinite number of forcing angular fre- 
quencies a,„ k in the primary. In general, the frequencies are 
different from zero so that the associated terms give rise to dy- 
namic tides. Terms in the tide-generating potential for which 
a,„ k = 0, on the other hand, give rise to static tides. We note 
that at least one such static tide exists for each spherical har- 
monic degree t. the tide generated by the term associated with 
k = m = 0. In binaries with spin-orbit resonances, additional 
static tides exist for non-zero values of k and m satisfying 
k/m = —fli/n. 

From the definition of the Fourier coefficients C(,,„± it fol- 
lows that the £ = 3 and £ = 4- terms in the expansion of the 
tide-generating potential contain additional factors Ri/a in 
comparison to the £ = 2 terms. The expansion of the tide- 
generating potential is therefore usually restricted to the dom- 
inant 1 = 2 terms. However, in the case of binaries in which 
the primary is close to or is filling its Roche lobe (e.g. inspi- 
ralling double white dwarfs and AM CVn binaries), the ratio 
Ri/a may become large enough so that the £ = 3 and £ = 4 
terms are no longer negligible. 

In the particular case of a binary with a circular orbit, the 
Fourier coefficients C£ ,„ ^ are different from zero only when 
k = -m. The non-zero coefficients are then given by 

_ (^-H)! pM. 



'(0) 



(5) 



^ For some I and m values the coefficients < increase to a local maxi- 
mum before decreasing with increasing values of k. 
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values of £, though in the apphcations we restrict ourselves to 
mass ratios ^ > 1 for which it is sufficient to only consider 
terms in the tide-generating potential associated with £ = 2. 

3. THE EQUATIONS GOVERNING FORCED 
OSCILLATIONS OF A SPHERICALLY SYMMETRIC 
STAR 

When the tidal force exerted by the companion is treated as 
a small perturbing force acting on a static, spherically sym- 
metric, equilibrium star, the equations governing the tidal dis- 
placement field ^tide and the perturbations of the star's mass 
density p, pressure P, and potential of self- gravity $ are ob- 
tained by perturbing and linearizing the equation of motion, 
the equation expressing the conservation of mass, the energy 
equation, and the Poisson equation. The perturbed and lin- 
earized equations take the form (e.g. Ledoux & Walraven 
1958) 



Fig. 1. — Absolute value of the non-zero Fourier coefficients C£ „, _„, as 
functions of tlie mass ratio q = M2/M1, for circular binaries with a Roche- 
lobe filling primary. Since C£,„_,„ = ci-,„ „,, only Fourier coefficients with 
positive values of m are shown. 

and the tide-generating potential reduces to 

GMi ^ ^ 

etide W (r, t) = -Etide — 2^ 2^ Ce,m-m 



1=2 , 



F;"(6i,0)exp[iOT(r2if-M)]. 



(6) 



Hence, for binaries with circular orbits, the expansion of the 
tide-generating potential becomes independent of time with 
respect to the co-rotating frame of reference when the star's 
rotation rate is synchronized with the orbital motion. The tide- 
generating potential then corresponds to the potential giving 
rise to equilibrium tides. The values of the coefficients q „, ^ 
for binaries with circular orbits are listed in Table [T] 

If the primary in a binary with a circular orbit fills its Roche 
lobe, its radius Ri is approximately equal to the volume- 
equivalent radius of its Roche lobe. In terms of the mass 
ratio q = M2/M1, the latter can be expressed as (Eggleton 
1983) 

^ ^ 0.49 g-^/^ 

a 0.6^-2/3 + ln (1+^-1/3)' ^' 

so that the Fourier coefficients ,„jt become a function of the 
mass ratio q. The variations of these coefficients as functions 
of q are shown in Fig. [T] For large mass ratios ^ 1, the 
second-degree Fourier coefficients C2.m.k are always at least an 
order of magnitude larger than the Fourier coefficients C2.,„.k 
and C4 so that the £ = 3 and £ = 4 terms can be neglected 
in the expansion of the tide-generating potential. However, 
when q ^1, the magnitude of the Fourier coefficients C3 i 
C3-1.1, and C4_o.o becomes comparable to that of the coeffi- 
cients C2,2.-2 and C2.-2,2- The other non-zero Fourier coef- 
ficients remain at least an order of magnitude smaller than 
C2,2.-2 and C2.-2,2- Since non-Roche-lobe filling primaries al- 
ways have Ri < the curves shown in Fig. [Upose upper 
limits on the Fourier coefficients „,_i in detached binaries 
with circular orbits. 

In the following sections, we derive a formalism to study 
quasi-static tides in close binaries. Given the potentially non- 
trivial role of the £ = 3 and £ = 4 terms, we allow for arbitrary 



tide 



Ptide 
P 

d_ 
dt 



= -V*tide- 

1 dp 



1 \ dP , 

-V^tide + ^^Pude 



p dr 

P P dr^' 



Pt'ide , 1 / 

p p dr 



(Ls-Dp 



d_Q 
dt 



tide 



V'*;Me=47rGp;i 



idei 



(8) 
(9) 

(10) 
(11) 



where 4* tide = ^'dde + Stide W is the total perturbation of the 
gravitational potential, the radial component of the tidal 
displacement field, dQ/dt the rate of change of thermal en- 
ergy (see Appendix |B] for details), and a prime on a quantity 
denotes the Eulerian perturbation of that quantity. The gener- 
alized isentropic coefficients Fi and are defined as 



and 



Fi = [(ainP)/(ainp)],, 



V^-\ = [(d\nT)/{d\np)]s, 



(12) 



(13) 



where S is the entropy, and T the temperature (e.g.. Cox & 
GiuH 1968). 

Following Savonije & Witte (2002), we furthermore ac- 
count for the effects of turbulent convection on the tides raised 
by the companion by adding a radial viscous force density 



/r.visc — -) r^ 



l_d_ 
pr~ dr 



d f 96,tide 



or \ 



dt 



(14) 



to the right-hand member of Eq. (O. In this equation, u is 
the coefficient of turbulent viscosity. The viscous force is dif- 
ferent from zero only in convective regions of the star. As 
outlined in Appendix|B] we do not incorporate effects of con- 
vection on the perturbation of the rate of change of thermal 
energy. 

Next, we separate the time and angular coordinates in 
Eqs. (l8]l-(fmi by expanding the tidal displacement field as 



£=2 I 



ek=- 



{r) d rii,„,,t{r) d 



86' 



:Y^\e,(f>) exp[i {a,„,tt-knT)] 



rsin6' dcf) 
(15) 
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Fourier coefficients cf ., 



TABLE 1 

i„ FOR BINARIES WITH CIRCULAR ORBITS. 



m = 


m = ±1 


m = ±2 


m = ±3 m = ±4 


1 = 2 -Ml 





1/8 




1 = 3 


(l/8)(Ri/fl) 





-(l/48)(i?,/fl) 


i = A (3/8)(Ri/«)2 





-(l/48)(RiA()2 


(l/348)(Ri/a)2 



and the perturbations of the stellar structure quantities as 

A t oo 
1=2 m=-^k=—oo 

X exp [i((T„,,jif-A;nr)] . (16) 

Substitution of the expansions in Eqs. (l8ll-(fTTli and replacing 
the 9 and ((> components of the equation of motion by the ra- 
dial component of the vorticity equation and the equation for 
the divergence of the transverse component of the tidal dis- 
placement field (see, e.g., Ledoux & Walraven 1958, Aizen- 
man & Smeyers 1977, Willems 2000) yields 



i.m.k 



.m,k ■ 



dr 

1 d 



dr p dr 



2 d(^(,„f: 

-^(Tm,k—-r Pr V—- 
pr~ dr \ dr 

P' 

^lm,k VLin^k = ^e.m,k + , 
P'l.m,k , dp 

1 y 'il.m,k - ~<^e.m,k, 

p p dr 



(17) 
(18) 
(19) 



^CFm.k 



— r n ^ t,e.m.k+^ I Oie.m.k 

P P dr 



]_d_ 

r^ dr 



_ (rw)p (dQ 
P \dt 

,d<i/ 



(20) 



i,m,k 



dr 



= 4^GpL,,. (21) 

Here, ae^m,k is the divergence of the tidal displacement field 
given by 



W . , . . £{£+1) 

i,k = — -r [r l;e.m.k ) o — Ve.nhk- 

r~ dr r- 



(22) 



In the next section, we derive approximate solutions to the 
above system of differential equations in the Umiting case of 
small forcing frequencies. For this purpose, it is convenient to 
use Eqs. ( fTSl l, ( fT$l ), and ( l20b to eliminate ^, p', and P' from 
Eq. OTJi. It follows that 

2 2 f ^ d-qi,m,k 



1 d 
pr^ dr 



+ io-„a--^— 1 pr V 
i 1 dp 



dr 

d£,e.m,k 
dr 



._.^<r,-i,(fV 

CT,„.* P dr \dt 



(23) 



i.m,k 



where A^^ = -g[{g/c^) + {dlnp/dr)] is the square of the Brunt- 
Vaisala frequency, g the gravity, and Cj = T\P/ p the square of 
the isentropic sound speed. 



Equations (fT7Ti-(l2ni must be supplemented with boundary 
conditions at the star's center and at the star's surface. At 
r = 0, we impose that the tidal displacement field remains fi- 
nite. At r = R\, we adopt zero-boundary conditions and im- 
pose the Lagrangian displacement of the pressure to vanish: 
(<5^')tide(^i) = 0. We furthermore impose the gravitational po- 
tential and its first derivative to be continuous at r = /?i , which 
is expressed by the condition 

d^emk\ ^+1 



dr 



= -etide(2^+ 1)^-2^ eg jn^k 

R, 



(24) 



(e.g. Polfliet & Smeyers 1990). Because of the non- 
homogeneous term in the right-hand member of this equation, 
the solutions to Eqs. (fT7]i-(l2Tli are proportional to the product 

^tideCe.w.k- 

A system of equations of the form of Eqs. (fT7l)- (l2T]) exists 
for each £, m, and k in the expansion of the tide-generating po- 
tential. The system of equations is complex due the presence 
of the convective damping term in the right-hand member of 
Eq. ( [TtI i and the perturbation of the rate of change of thermal 
energy in the right-hand member of Eq. ( |20] |. The solutions 
will therefore show a phase shift with respect to those found 
in the adiabatic approximation. We furthermore note that the 
solution to Eqs. (fT7]i-(l2ni associated with the forcing angular 
frequency -cTm /t is the complex conjugate of the solution as- 
sociated with the forcing angular frequency (t„, The two so- 
lutions therefore have the same amplitude, but opposite phase 
shifts. 

4. THE QUASI-STATIC TIDE APPROXIMATION 

Smeyers (1997) and Smeyers & Willems (1998) derived ap- 
proximate solutions to Eqs. (fT7b- (l2T]) for low-frequency dy- 
namic tides valid to 0(0-^,^). The authors adopted the adi- 
abatic approximation and found the radial component of the 
tidal displacement field to be described by a non-oscillatory 
term of Oia'^^^i,) and an oscillatory term of 0{(Jj^i^) (see also 
Zahn 1975). No term of O(am.k) occurs in the adiabatic ap- 
proximation. Because of the oscillatory nature of the solu- 
tions at order Oiaj^^ i^), a multi-variable perturbation method 
was adopted to account for the large second derivatives asso- 
ciated with the rapid variations of the radial component of the 
tidal displacement as a function of the radial coordinate r. 

When convective damping and the perturbation of the rate 
of change of thermal energy are taken into account, terms of 
0{a,„ k) appear in the low-frequency approximations of the 
solutions to Eqs. (fT7b- (l2T]) . We therefore look for solutions 
valid to 0{am,k) in the limiting case where the forcing angu- 
lar frequency is small enough to treat the tides as quasi-static 
in the frame of reference co-rotating with the star In this 
limiting case, the solutions do not show rapid variations in the 
radial direction, so that no multi-variable perturbation method 
is required to derive the approximate solutions. Our approach 
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TABLE 2 
Units of physical quantities 



Quantity 


Unit 


Quantity 


Unit 


f 




r 








V 




P 


Mi/{4nR]) 


P 


GM\I {A-kR\) 


$ 


GMi/Ri 


etidc W 


GMi/Ri 


V 


(GMiR|)l/2 


dQidt 





is similar in nature as that of Campbell (1984), except that 
we incorporate the effects of convective damping and allow 
for arbitrary orbital eccentricities. We also make the deriva- 
tion more transparent by using the total perturbation of the 
gravitational potential (which is of 0(crj', j)) rather than the 
tidal velocity field (which is of 0((Tm,A)) as the main depen- 
dent variable. 

To derive approximate solutions to the system of equations 
governing tides in close binary components in the limiting 
case where \a,„,k\ <C (GMi//?j)'/^, it is convenient to pass on 
to dimensionless quantities by expressing the physical quan- 
tities in the units listed in Table |2l After passing on to dimen- 
sionless quantities, Eqs. ( l20l i and ( l23T i can be written in the 
form 



:,A- = 



1 d 



drii,, 



pr 
C I dp 



pp' V 



dr 

d£,e.m,k 
dr 



+1 -^(r3- 

crm,k P dr 



' \ dt 



(25) 



l.m,k 



and 



P',„k^dP 

— r n +i I OLijn^k 

P P dr 

C (Ti-\)pfdQ\ 
dt 



(26) 



e..m,k 



In these equations, C is the ratio of the star's dynamic time 
scale to its Helmholtz-Kelvin time scale: 



1/2 



/ 



( GM\ 
\RiLi 



(27) 



In general, C <C cr,„_jt so that the terms containing the Eulerian 
perturbation of the rate of change of thermal energy in the 
right-hand member of Eqs. dZSl l and ( |26] | are non-negligible 
only near the star's surface where the factors {l/p){dp/dr) 
and p/P become large. We therefore assume these terms to 
be of 0(a„,^k) at least in some region near the star's surface 
and, following Willems et al. (2003), set 

C 

'^m,k 



C'e. 



(28) 



Next, we expand the components of the tidal displacement 
field and the perturbed stellar structure quantities in series of 
the form 



fiM.k(r) = f!'Ur) + a„,,n:>Jr) + 



<i) 



(29) 



For brevity, we omit the subscripts £, m, and k from the com- 
ponents of the tidal displacement field and the perturbed stel- 
lar structure quantities for the remainder of this section. 



After substitution of the expansions for the components of 
the tidal displacement field and the perturbed stellar structure 
quantities, Eq. dZST l, at ©(cr"), yields 



.(0) : 



:0 



(30) 



in regions of the star where A^^ 0. From Eqs. ( fTSl ), ( fT9b , 
and (l26b . it then follows that 



^(0) = . 



pm 
p 



dp MO) 



p'(0) = _^^(0) 

dr 



= -- — . 



(31) 



Substituting these equations into Eq. (|2TI) and making 
use of the equation of hydrostatic equilibrium, we derive a 
second-order differential equation for given by 



^2^(0) 2 £/*(0) 

■ + - 



dr^ 



dr 



I dp £(i+l) 
g dr r^ 



= 0. (32) 



This equation is equivalent to the equation of Clairaut for 
the radial component of the tidal displacement field usually 
derived in the framework of the theory of equilibrium tides 
(Sterne 1939). At C'(cr''), the boundary condition given by 
Eq. ( |24] | takes the form 



) + f ^+1-^ ) =-etide(2£+l)Q,,„i, (33) 



d^' 
dr 

where ps and gs are the mass density and gravity at the sur- 
face of the unperturbed star, respectively. The solution to 
Eq. ( |32] | that remains finite at the star's center and satisfies 
the surface boundary condition is real and proportional to the 

product £tideQ,ml- 

At 0{a), Eq. i25[ for the divergence of the tidal displace- 
ment field yields 



g 1 d 
c^N^ pr^ dr 



d^' 
pr V 



(0) 



+ iC' (Fs-l) 



dr 

g 1 dp (dQ 



clN^ p dr \ dt 



/(O) 



(34) 



in regions of the star where A^- ^ 0. After substitution of this 
solution into Eqs. ( fTST i, ( fT9] ), and ( l26b , we find expressions for 
the total perturbation of the gravitational potential, the Eule- 
rian perturbations of the mass density and pressure, and the 
radial component of the tidal displacement field given by 



P'{i) 



^/(I) __^^(I)_ 



p-- ' = - 



dr 



pa 



(1) 



y(i) . 



^.(i) 
dr^ 



iC (r3-i)p 



dt ) 



/(()) 



(35) 



(36) 



(37) 



e(" = — +i 



-r\prv- 



dC 



(0) 



g pr^N^ dr 



dr 



'(()) 



(38) 
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Finally, proceeding in a similar way as for the derivation 
of Eq. (|32] |. leads to a second-order differential equation for 



■ + - ■ 



1 



dr 
d 



gr~ dr\^ dr 



I dp £(£+l) 
g dr r^ 



(39) 



The associated boundary condition expressing the continuity 
of the gravitational potential and its first derivative takes the 
form 



dr 



1 



N^r^ dr 



+ iC' 









-) 

pr 




P8 


fdQ 



c^N^ V dt 



'(0) 



(40) 



At 0{(j), the differential equation and surface boundary con- 
dition for the total perturbation of the gravitational potential 
are both complex. The solution therefore shows a phase shift 
with respect to the tide-generating potential, and thus also 
with respect to the instantaneous position of the companion. 
The total perturbation of the gravitational potential further- 
more depends on the perturbation of the rate of change of 
thermal energy only through the boundary condition express- 
ing the continuity of the gravitational potential and its first 
derivative. In AppendixiBl we show that the lowest-order per- 
turbation of the rate of change of thermal energy, {dQ/ dt)'^^^\ 
is determined by the lowest-order solution for the total pertur- 
bation of the gravitational potential, 'I'^^', and its first deriva- 
tive, d'^^^^'/dr. The solutions at 0{<j) are therefore also pro- 
portional to the product etideC£,m.A. A semi-analytical solution 
method for the system of differential equations composed of 
Eqs. ( l32b and ( l39b and their associated boundary conditions 
is outlined in Appendix Icl From the solution, it follows that 
^'^''(r) is purely imaginary. 

5. ORBITAL EVOLUTION 

The tidal distortion of a star perturbs the spherical symme- 
try of it's external gravitational field, which in turn perturbs 
the motion of the companion from a pure Keplerian orbit. We 
study the perturbed motion in the framework of the theory of 
osculating elements in celestial mechanics (e.g. Sterne 1960, 
Brouwer & Clemence 1961, Fitzpatrick 1970). In this frame- 
work, the rates of change of the orbital semi-major axis and 
eccentricity due a star's tidal distortion are given by 

da 2 dTZ 
dt n^a dr ' 



de 
'dt 



1 



n dr 



+ (1- 



2N1/2 dn 



(41) 
(42) 



were XS is the longitude of the periastron, and 7?. is a perturb- 
ing function related to the Eulerian perturbation of the star's 
external gravitational potential $',('% t) as 

Mx+Mi 



TZ{u,v,t) ■■ 



'2' 



-Qit,t 



(43) 



(Smeyers et al. 1991). The Eulerian perturbation of the ex- 
ternal gravitational potential is evaluated at the position of the 



companion, r = u, 9 = tt/2, and cf) = v- flit, where u is the 
distance between the stars. 

Since 7^ is a function of the distance u and the true anomaly 
V, we transform the partial derivatives with respect to t and (53 
in Eqs. dTO and (l42l i into partial derivatives with respect to u 
and V. By the use of the equations (e.g. Fitzpatrick 1970) 

du nae ") 

smv, 



dr 
dv 



1- 



n 



(1- 



.3/2 



(1 +e COS v)" 



(44) 



the rates of change of the orbital semi-major a and eccentricity 
e due to the tidal distortion of a binary component then take 
the form 



da 

dt nail 



de 
It 



2e 



,1/2 



dn {l+ecosvf dn 
asmv — + — J- rr— 



(45) 



na^ 

{ . dn I 

X < fl sm V — — I- - 
du e 



(1 +e cos v)" 

r^^2 



1 



) ■ (46) 



The Eulerian perturbation of the external gravitational po- 
tential due to the primary's tidal distortion is a solution of the 
equation of Laplace. The solution which tends to zero at in- 
finity can be cast in the form 



(47) 



l=Q m=-tk=-oo 

xF;'(0,(/))exp[i(cr,„,if-^«r)] , 

where the A^.mjt are constants determined by the condition 
that the Eulerian perturbation of the gravitational potential be 
continuous at the star's surface. 

From the definition of the total perturbation of the gravi- 
tational potential, ^Ptye = ^t'jde + £tide^> and Expansions ([3j 
and ( fT6l ), it follows that the Eulerian perturbation of the star's 
potential of self-gravity can be decomposed in terms of spher- 
ical harmonics and Fourier series as 



'i>:ide(^o=EEE 



.ml(r) + Etide — Q,m.A- 



1=2 m=-lk=-c 

xy;'(6l,</.)exp[i(cr„af-/t«T)] . (48) 

Continuity of the gravitational potential at the star's surface 
thus requires 



A 
with 



GMi 

k = -^tide — — Cij 



Ri 



Ft, 



for ^ = 2, 3,4, 
otherwise. 



l.m.k' 



(Ri) 



+ 1 



(49) 



(50) 



The quantities Fi,,n.k are dimensionless and measure the re- 
sponse of the star to the various forcing angular frequencies 
(j„, J appearing in the expansion of the tide-generating poten- 
tial. Since 5'£,m,A- oc EtideQ.m.i^ the quantities are independent 
of the product etide ce^„,,k- In the limiting case of long orbital 
and rotational periods, the forcing angular frequencies (t„,i 
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all tend to zero and the constants F(,„,± tend to the classical 
apsidal-motion constants ke determined in the framework of 
the theory of static tides (Smeyers & Willems 2001). 

For tides with non-zero forcing frequencies, the quantities 
Fi.m.k are complex. It is therefore convenient to write them in 
polar form as 



Fi,m,k = \Fij„^k\ exp [i^e.m.k] 



(51) 



Since the solution to Eqs. (fT7]i-(l2ni associated with the forc- 
ing angular frequency -(J,„,k has the same amplitude, but op- 
posite phase as the solution associated with the forcing an- 
gular frequency (J„,,k, the magnitude and phase of the quanti- 
ties Fe^„,^k obey the properties |f>.-,„.-i| = |fi,m,/i|, and ji-m-k = 
-li.m.k- By means of these properties and the symmetry prop- 
erties of the coefficients c^ ,,, i. Expansion WT\ for the Eule- 
rian perturbation of the external gravitational potential can be 
written in real form by combining the terms associated with 
the forcing angular frequency a-m-k with the terms associated 
with the forcing angular frequency amy. 



TZ{u,v,t) = 4- 



G{Mi+M2) M: 



M2 \ ^ \ ^ ^— \ / Ri 

M, 2^ 2^ 2^ \ n 



M 

-(£+1) 



1=2 m=-e k=0 



PPiO) 



y-Hi^m,kCi,m,k\F£,m,k\[-) COS [mv + kM + '-fi^„j,k) , (52) 

where 



(53) 



Finally, by transforming the time derivatives in Eqs. ( |45] l 
and ( l46b into derivatives with respect to the mean anomaly 
and averaging over one revolution of the companion, we de- 
rive the equations for the rates of secular change of the orbital 
semi-major axis and eccentricity due to the tidal distortion of 
a binary component to be 





= 1/2, 








Ke,m,0 


= 


for 


-l<m 


<-l, 




= 1 


for 


I <m 




l^l,m.k 


= 1 


for 


-e<m 


<£nndk>l 



Stt M2 



'EEE 



e=2 , 



X Kf jia \Ftj,a\ sin7f,„,,(tG^„,_^(e) 



-e k=o 

<2) 



e+3 



+00 / „ X e+3 
R\ 



(54) 



_ Stt M2 

X K.(j„,k \F(j,a \ sin7f,„,,A-G^'^i_<.(e) 



1=2 m=-e k=0 
-(3) 



(55) 



where 



Gflkie)-- 



1 



Ctm.kPi"\0) 



X ■ 



TT 



(£+l)e I (l+ecosv) sin(mv-l-fcM) sinvc/v 




-m I (1+e cosv)^^^ cosimv+kM)dv 
'0 



(56) 
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Fig. 2. — Absolute value of the non-zero coefficients G?' (e) as func- 
tions of the mass ratio q — M^jMi, for circular binaries with a Roche-lobe 



filhng primary. Since G^^^ = - 
positive values of m are shown. 



-a) 

■'1,-m.m 



(e), only the coefficients with 



G?L(^) = 



rQ.„,tri"''(0) 



X— ■^(£-l-l)e / (1-i-ecosv) sin(mv-l-^M)sinviiv 
-m I (1-l-ecosv)''"' [(1 -l-e cos v)"- (l -e^)] 



X cos(m V + kM) dv 



(57) 



The coefficients Gfl^/^ie) and Gflji,(e) are the same as those 
derived by Willems et al. (2003) for the rates of change of 
the orbital semi-major axis and eccentricity due to resonances 
between dynamic tides and free oscillation modes of close 
binary components. They are functions of the orbital eccen- 
tricity e and, through the Fourier coefficients cgjn.k, are pro- 
portional to the ratio {Ri/aY~^. They obey the properties of 



asymetry G^l„,_,(.) = -GZ,ki^) and G-i,„,_,(.) = -G-;,(e), 
and are different from zero only when the corresponding co- 
efficients ci,m.k are different from zero (see Sect.©. The coef- 
ficients G^^Q g(e) and G^^'g o(e) are furthermore identically zero 

for all orbital eccentricities e. Similar coefficients G£'j,^(e) 
exist for the rates of secular change of the position of the pe- 
riastron (see Smeyers et al. 1998, Willems 2000, Willems et 
al. 2003). 

In the particular case of a binary with a circular orbit, the 
coefficients G^^\^ ^ are different from zero only when k = —m 
and m ^Q. The non-zero coefficients take the values 



-(2) 



-(3) 



-(3) 



G 



!(2) 



(58) 



and are listed in Table [3] If the primary in a binary with a 
circular orbit furthermore fills its Roche lobe, the coefficients 
k become functions of the binary mass ratio q. The vari- 
ations of these coefficients as a function of the binary mass 
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Coefficients Gf]^ 


TABLE 3 










„j FOR BINARIES WITH CIRCULAR ORBITS. 






m = 


m = ± 1 


m = ±2 


m = ±3 


m = ±4 


1 = 2 








T3/2 






i = -i 





T(3/8)(Ri/fl) 





T(15/8)(filA') 




e = 4 








T(15/24)(R|/a)2 





T(35/16)(Ri/a)2 



ratio are shown in Fig. [2] For mass ratios ^ ^ 1, the second- a magnitude 



(2) 

degree coefficients ^ are always an order of magnitude 
or more larger than the coefficients Cfl^ and Cfl^ so that 
the £=3 and £ = 4 terms can be neglected in the expansion 
for the rate of secular change of the orbital semi-major axis. 

(2) 
3.m.k 



For mass ratios < 1, the magnitude of the coefficients G 



and GJj^*, j can become comparable to the magnitude of the 



coefficients G2 I, f.- Since non-Roche-lobe filling primaries al- 
ways have Ri < Rl,i, the curves shown in Fig. |2]pose upper 
limits on the coefficients G?i in detached binaries with cir- 



'i,m,k 
(3) 

cular orbits. The coefficients Gy^j^ are all identically zero for 
binaries with circular orbits (Willems 2000). 

In the limiting case where all relevant forcing angular fre- 
quencies a,„_k are small compared to the inverse of the star's 
dynamical time scale, the quantities /v_,„jt can be determined 
by means of the quasi-static tide solutions derived in §|4l 



1 

Fe,m,k = -- 



(1) 



GMi 



■ + 1 



where ^f,*„i(^i) is real and ^^,',,^(^1) is imaginary. The 
phase angles 7^ ,,„,Ar are then determined by 



Slide Q,ml- 



(I) 



(59) 



tan ji^,n,k = ■ 



Slide C£ in.k 



^1 



(0) 
e.m.k 



Slide Ci,m.k 



(Ri) GMi 
+ ■ 



(60) 



When dissipative effects are small, tan7£,„,jt « "fi,m,k, so that 
the phase angles 7^ m.i are proportional to the forcing angular 
frequencies am,k- This proportionality is usually assumed in 
the context of the weak friction approximation in tidal evolu- 
tion theory (see, e.g., Alexander 1973, Hut 1981, Ruymaek- 
ers 1992). In Appendix iDl we show that in the limiting case 
of small forcing frequencies and weak damping, Eqs. ( |54] | 
and ( |55] l are equivalent to the equations derived by, e.g., Zahn 
(1977, 1978), Hut (1981), and Ruymaekers (1992). 

6. STELLAR SPIN EVOLUTION 

Due to the phase shift between the perturbation of the grav- 
itational potential and the position of the companion induced 
by tidal energy dissipation, the companion exerts a torque T 
on the tidally distorted star. The torque is determined from 
Newton's law of action and reaction as the opposite of the 
torque exerted by the tidally distorted star on the companion 
due to the perturbation of the star's external gravitational po- 
tential: 



f = M2(rx V$^), 



(61) 



where is to be evaluated at the position of the companion. 
The tidal torque is perpendicular to the orbital plane and has 



T: 



87r 

X mP\ 
u 
a 



gm\mI 

Ml +M2 



1/2 4 f 00 . ^ ^+3 

— fli/^V W ( ill 
Ml p 

' e=2 m=-l k=0 ^ 



' iO)Ke.m,kCe,m.k \Fe,m,k\ 

sin (mv+kM+jij, 



(62) 



Since we are interested in the long-term secular effects of 
the tidal torque on the tidally distorted star, we average the 
torque over one revolution of the companion. It follows that 



with 



8tt f GM\M: 



1/2 



Ml 



I sin7f.,„_A:G^'*^j.(e) 



EEE 



1=2 m=-e k=0 



e+3 



.(4) ii+\m\y. 



{£-\m\)\ 



-d-i) 



(63) 



(64) 



The coefficients G^f'^ f.(e) are different from zero only for non- 
axisymmetric (m ^ 0) tides and have the same sign as the 
azimuthal number m. They are related to the coefficients 
G?' (e)and Gf„.{e) as 



(1 



.^2)1/2 



le 



■G%Ae) 



(65) 



so that their properties can be derived either from the prop- 
erties of the coefficients c^„,_jt or from the properties of the 

(3) c„\ 'pjjg values of the coeffi- 



coefficients G^'^\^ ^ 

-.(4) 



(^) and G^X ,,(.). 



cients G\ „, ^ for binaries with circular orbits are listed in Ta- 



ble|4] The variations of the coefficients G'"'^\ for binaries with 
circular orbits and Roche-lobe filling primaries are shown in 
Fig. [3] as functions of the binary mass ratio. 

The tidal torque T exerted by the companion affects the 
rotational angular velocity of the tidally distorted star As- 
suming the star rotates as a rigid body and neglecting any 
perturbations of the star's moment of inertia due to the tidal 
distortion, the rate of change of the rotational angular velocity 
ill is related to the tidal torque as 



d^i ^ 



(66) 



where li is the star's moment of inertia with respect to its 
rotation axis. Consequently, the rate of secular change of the 
rotational angular velocity ili is given by 

GM\M\\"^ M2a^''^ 

Ml h 2^ 2^ 2^ 



dVLi 
dt 



87r 



Ml +M2 



=2 m=-l k=0 



Hi,m,k \Fi^m,k\ sm^i_„aG[l^(e) 



(67) 
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TABLE 4 

Coefficients cfl for binaries with circular orbits. 



m = 



m = ±l 



m = ±2 



m = ±3 



m = ±4 



1 = 2 

^ = ^ 

^ = 4 



±3/4 
±(3/16)(Ri/a) 

±(5/16)(Ri/a)^ 



±(15/16)(Ri/fl) 

±(35/32)(fii/a)^ 



II 
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j'/' (e) as func- 
for circular binaries with a Roche-lobe 
filling primary. Since — ~^'"T-m nS^"^' ""'^ coefficients with 

positive values of m are shown. 



Fig. 3. — Absolute value of the non-zero coefficients G\ 
tions of the mass ratio q = Mi/ Mi 

'i 



7. ASTROPHYSICAL RELEVANCE 

We apply the formaHsm presented in the previous sections 
to a Q.3Mq hehum white dwarf model representative of an 
isolated white dwarf with a radius of O.O18/?0 and an effec- 
tive temperature of 3590 K. Other relevant model properties 
are summarized in Table |5] and details on the model input 
physics can be found in Deloye et al. (2007). We particu- 
larly note that the model has a thin convection zone near the 
stellar surface characterized by a frequency-dependent turbu- 
lent viscosity coefficient 



1 + 



In 



(68) 



Here, L is the mixing length, Tconv = \N^\~^^^ the convective 
turnover time scale (Terquem et al. 1998, Savonije & Witte 
2002), and s an integer constant. The mixing length is as- 
sumed to be twice the local pressure scale height. 

The factor between square brackets in Eq. ( l68l l is a reduc- 
tion factor to account for the decreased efficiency of convec- 
tive damping when the tidal period P„,,k = 2TT/a„,,k is shorter 
than the convective turnover time scale Tconv Because of 
the dependence of the reduction factor on the tidal forcing 
frequency (T,„_jt, the turbulent viscosity coefficient is differ- 
ent for tides generated by different terms in Expansion (|3]l 
for the tide-generating potential (see also Zahn 2008). Zahn 
(1966) proposed a reduction of the turbulent viscosity coeffi- 
cient characterized by s=\, while Goldreich & Keeley (1977) 
argued for s = 2 (see also Goldman & Mazeh 1991, Goodman 
& Oh 1997). More recently, Penev et al. (2007) calculated 
the effective turbulent viscosity as a function of the tidal forc- 



TABLE 5 

White dwarf model properties. 



Quantity 



Value 



Ml 

U 

b 

Tdyn.l 

l/fHK,l 
Chemical Composition 



O.3M0 
O.O18i?0 
3590 K 
4.7 X IQ-^ Lq 
6.9 s 

3.4 X lO'V 

6.5 X 10-20 
helium 



' Effective temperature'' Dynamic time scale'^ Helmholtz-Kelvin time scale 



6x10' 
5x10' 
V 4x10' 
^ 3x10' 

E 

2x10' 
1x10' 







; CT„.k=o.oi/T,,„,, 




T ct„.,=0.1/t,,„,, 




/ ' 




i / ' 





0.996 



0.997 



0.998 
r/R, 



0.999 



1.000 



Fig. 4. — Variation of the turbulent viscosity coefficient Vi„ t in the surface 
convection zone of a O.3M0 hehum white dwarf model for different forcing 
angular frequencies cr,„ i^. The a„, t = line represents the maximum attain- 
able turbulent viscosity. The dashed and dotted lines represent the reduced 
turbulent viscosity coefficient for (j,„ ■ - n ni -^-1 
suming s=\. 



ing frequency using 3-D numerical simulations and found a 
reduction factor that closely matched the s = \ prescription 
proposed by Zahn (1966). In our calculations, we therefore 
adopt s=\. Test calculations with s = 2 show that the choice 
of s affects our results by less than an order of magnitude. 

In Fig. m we show the variations of the turbulent viscos- 
ity coefficient v,„^k in the O.3M0 helium white dwarf model 
as a function of the normalized radial coordinate r/R\. The 
maximum attainable turbulent viscosity is indicated by the 
Cm.i: = line. The dashed and dotted lines show the re- 
duced turbulent viscosity for s= \ and tidal forcing frequen- 
cies cr„,,t = 0.01rjy'„ , and (T,„,t = O.lr^^'^ j, respectively The 
reduction of the turbulent viscosity coefficient is most promi- 
nent near the base of the convection zone where the convective 
turnover time scale is longest. 

For non-degenerate stars dissipation of tidal energy through 
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0.0 0.2 0.4 0.6 
r/R; 



0.8 0.995 



1.000 



r/R, 



Fig. 5. — Order of magnitude of the convective and radiative damping 
terms in the system of differential equations governing quasi-static tides in a 
0.3 Mq helium white dwarf model. The terms are calculated for £ = 2 and 
are rendered dimensionless by expressing the physical quantities in the units 
listed in table |2] The convective damping term is calculated for a forcing 
angular frequency o;„ ii = 0. 

quasi-Static tides is dominated by convective damping (Zahn 
1977). This still holds true in white dwarfs, as illustrated in 
Fig. |5] where the variations of the radiative and convective 
damping terms 



1 d 
— — [ pr-iy„a- 
ar 



(0) 



dr 



and C 



dQ 
dt 



are shown as functions of the normalized radial coordinate 
r/Ri. The terms are obtained by numerically integrating dif- 
ferential Eq. ( |32] | for £ = 2, and using Eqs. (l3Tl i and dBlSI l to 
determine £,^"\r) and {dQ/dt)^°^ from the lowest-order total 
perturbation of the gravitational potential ^/("'(r). For the de- 
termination of the turbulent viscosity term i',,, /;, the forcing 
angular frequency (T„,,k was set equal to zero. The terms are 
furthermore rendered dimensionless by expressing the physi- 
cal quantities in the units listed in Table |2] For convenience, 
we also divided the terms by the scaling factor etideCc,,,,.*, so 
that the curves shown in the figure are independent of the bi- 
nary orbital period, eccentricity, and companion mass. In the 
convection zone, the turbulent damping term exceeds the ra- 
diative damping term by more than 10 orders of magnitude. 
We could therefore safely have set the constant C' equal to 
zero in the derivation of the quasi-static tide solutions pre- 
sented in §|4] 

In the following subsections we calculate the tidal distortion 
and the orbital evolution time scales due to quasi-static tides 
for both detached and semi-detached white dwarf binaries and 
for different ranges of binary component masses, orbital peri- 
ods, orbital eccentricities, and white dwarf rotation rates. We 
only consider mass ratios q> 1, so that £=3 and £ = 4 terms 
can be neglected in the expansion of the tide-generating po- 
tential (see §|2ll. 

7.1. Detached binaries 

We first turn our attention to tidal interactions in detached 
white dwarf binaries. Given the strong dependence of tidal 
effects on the ratio of the white dwarf radius to the orbital 
semi-major axis, we focus on short-period binaries with or- 
bital frequencies relevant to the Laser Interferometer Space 



Antenna, LISA (10~'*-10~' Hz). In all cases, the companion 
star is assumed to be more compact than the white dwarf so 
that we do not have to worry about Roche-lobe overflow from 
the companion. As indicated above, we also restrict ourselves 
to the dominant quadrupole tides generated by the £ = 2 terms 
in Expansion (O of the tide-generating potential. 

To calculate the tidal distortion of the 0.3Mq helium white 
dwarf model, we solve Eqs. (|32] | and ( [39] l for ^'^'^'(r) and 
^''"(r) with their respective boundary conditions for each 
non-zero and non-negligible term in Expansion (O of the tide- 
generating potential. Equations ( |32] | and ( |39] l are solved using 
a variable step Runge-Kutta integrator, following the semi- 
analytical procedure outlined in Appendix|C] The radial com- 
ponent of the tidal displacement field and the Eulerian pertur- 
bations of the stellar structure quantities are then determined 
from Eqs. (EB, (p3-(l38]l, and (|A6ll-(|A7ll- We recall that the 
inclusion of dissipative effects renders the solutions complex 
and introduces a phase shift between the tidal perturbations 
and the tide-generating potential. The amplitude of the tidal 
perturbations is dominated by the real part of the solutions. 

The total tidal displacement field and perturbations of the 
stellar structure quantities are obtained by adding the solu- 
tions associated with all non-zero and non-negligible terms 
in the expansion of the tide-generating potential. For in- 
stance, the total perturbation of the gravitational potential due 
to quasi-static tides is given by 

4 e oo 

*t.de(r,f) = 2^ ^ ^ Kf,„a |*£.„a('-)|^]'"'(cos0) 

e=2 m=-l k=0 

X COS [m0+(Tm,<:f-A:nr+T*^___j,(r)] , (69) 

where 



T~*f.„,.i('') = arctan 



(70) 



The total radial component of the tidal displacement field and 
the total Eulerian perturbations of the stellar structure quanti- 
ties are obtained similarly. 

In Fig. |6l the amplitudes of the radial component of the 
tidal displacement field and the Eulerian perturbations of the 
mass density, pressure, and temperature are shown as func- 
tions of the normalized radial coordinate r/Ri. The Q.3Mq 
helium white dwarf is assumed to have an equal mass com- 
panion (^ = 1) in a circular binary with orbital periods rang- 
ing from 3 to 60 minutes. The white dwarf rotation period 
is assumed to be 100 hours. As expected, the amplitudes of 
the perturbations increase outward and are larger for shorter 
orbital periods. The radial component of the tidal displace- 
ment field at the star's surface in particular decreases from 
2 X lQ~^Ri to 4 X 10~^7?i when the orbital period increases 
from 3 to 60 min. 

Next, we consider the tidal evolution time scales f„ = \d/a\ 
and r[2, = |r2i/f7i| for the orbital semi-major axis and the 
white dwarf rotational angular velocity. The time scales are 
shown in Fig.|7]as a function of the orbital period for a 0.3Mo 
helium white dwarf in a circular orbit around a O.3M0 or 
10^ Mq point-mass companion^. The different curves corre- 
spond to different white dwarf rotation periods ranging from 
1 to lOOOhr. 

' We note that at the shortest orbital periods shown (Po,i, ^ 3 min), the 
forcing angular frequencies cr^ ii can be close to O.Srdyn.i, stretching the ap- 
plicability of the apphed perturbation theory. 
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Fig. 6. — Amplitude of the radial component ^ of the tidal displacement 
field and the Eulerian perturbations of the mass density p, pressure P, and 
temperature T for quasi-static tides associated with 1 = 2 ma 0.3 Mq helium 
white dwarf model. The white dwarf is assumed to have a rotation period of 
100 hr and to have a O.SMq point-mass companion in a circular orbit. From 
top to bottom, the different curves in each panel correspond to orbital periods 
= 3, 6, 13, 28, and 60 min. 



Fig. 7. — Tidal evolution time scales for the rates of secular change of the 
orbital semi-major axis a and the white dwarf rotational angular velocity Ci 
for a O.3M0 helium white dwarf in a circular binary. Solid and dotted lines 
represent time scales for point-mass companions of O.3M0 and 10^ Mq, re- 
spectively. From top to bottom, the solid and dotted lines in each panel cor- 
respond to white dwarf rotation periods of 1, 10, 100, and lOOOhr. In the 
top panel, the curves associated with rotation periods longer than lOhr are 
indistinguishable. The grey horizontal fine represents the age of an isolated 
0.3Mq helium white dwarf when it has cooled to a temperature of 3590 K. 




Fig. 8. — Time scales for the rate of secular change of the orbital semi- 
major axis a due to gravitational wave emission for a 0.3Mq white dwarf 
with aO.3M0 or lO'M© companion in a circular orbit. 



The time scales for the rate of secular change of the orbital 
semi-major axis are longer than a Hubble time for all binary 
configurations considered. They increase with increasing bi- 
nary companion mass due to the associated increase of the 
orbital semi-major axis for a given orbital period. For compar- 
ison, the time scales of orbital evolution due to gravitational 
wave emission determined from the Peters (1964) equations 
are shown in Fig. [8] for the same binary component masses 
and orbital period range as used in Fig. |2l The time scales 
of orbital evolution due to convective damping of quasi-static 
tides are several orders of magnitude larger than those due to 
gravitational radiation, so that the effects of quasi-static tides 
in white dwarfs can safely be neglected when predicting grav- 



itational wave signals from white dwarf binaries. 

The time scales for the rate of secular change of the white 
dwarf rotational angular velocity are at least two orders of 
magnitude shorter than the time scales for the rates of sec- 
ular change of the orbital semi-major axis. For a given or- 
bital period, they decrease significantly with increasing white 
dwarf rotation period or, equivalently, with increasing degree 
of asynchronism. For white dwarf rotation periods of lOOOhr, 
the time scales become shorter than the age of an isolated 
O.3M0 helium white dwarf of 3590 K for orbital periods be- 
low ~ 20-25 min, depending on the white dwarf companion 
mass. However, for all binary configurations considered, the 



12 



Willems, Deloye, & Kalogera 




10 20 30 40 



Fig. 9. — Tidal evolution time scales for the rates of secular change of 
the orbital semi-major axis a, the orbital eccentricity e, and the white dwarf 
rotational angular velocity Ci for a 0.3Mq helium white dwarf in a binary 
with orbital eccentricity e = 0.3. The different hnes have the same meaning 
as in Fig. |7] In the top two panels, the lines associated with rotation periods 
longer than lOhr are indistinguishable. 

white dwarf spin-up time scales are still considerably longer 
than the orbital evolution time scales due to gravitational radi- 
ation. Quasi-static tides will therefore not be able to spin the 
white dwarf up fast enough to reach a synchronous rotation 
rate. 

The effects of the orbital eccentricity on the tidal evolution 
time scales are illustrated in Fig. |9] where the tidal evolution 
time scales = \d/a\, = \e/e\, and f^, = | are shown 

for a 0.3 Mq helium white dwarf orbiting a 0.3Mq or 10^ Mq 
point-mass companion with an orbital eccentricity e = 0.3. For 
a given orbital period, the time scales are shorter than those 
for a circular binary due to the stronger tidal interactions tak- 
ing place at the periastron of the binary orbit. However, the 
time scales for the rates of secular change of the orbital semi- 
major axis and eccentricity remain longer than a Hubble time 
for all binary configurations considered. Depending on the 
binary companion mass, the time scales for the rate of sec- 
ular change of the white dwarf's rotational angular velocity 
can become shorter than the age of an isolated 0.3 Mq he- 
lium white dwarf of 3590 K for orbital periods Porb ^ 30 min, 
provided the initial degree of asynchronism at periastron is 
sufficiently high. 

7.2. Mass-transferring binaries 




10° 10' 10^ 10^ 10" 10^ 10^ 



Fig. 10. — Tidal evolution time scales for the rates of secular change of 
the orbital semi-major axis a and the white dwarf rotational angular velocity 
C 1 as a function of the binary mass ratio q = M2 /M[, for a Roche-lobe fining 
O.3M0 helium white dwarf in a circular binary. From top to bottom, the lines 
in the bottom panel correspond to white dwaif rotation periods of 1, 10, 100, 
and lOOOhr. In the top panel, lines associated with the rotation periods are 
indistinguishable. The grey horizontal hne represents the age of an isolated 
0.3 Mq helium white dwarf when it has cooled to a temperature of 3590 K. 

Next, we consider tidal interactions for binaries in which 
the considered O.3M0 helium white dwarf fills its Roche lobe. 
We ignore any coupling between mass transfer and tides and 
assume the radius of the white dwarf to be exactly equal to 
the radius of its Roche lobe. The orbital separation is then 
fully determined by the mass ratio by means of Eq. ^ for the 
volume-equivalent radius of the white dwarf's Roche lobe. 
We refrain from comparing the tidal evolution time scales 
with the time scales of orbital evolution due to mass transfer 
since the determination of the latter requires detailed track- 
ing of the response of the white dwarf to mass loss (see, e.g., 
Deloye & Taam 2006, Deloye et al. 2007). 

The time scales tg = \d/a\ and fji, = |r2i/ili| for the rates 
of secular change of the orbital semi-major axis and the white 
dwarf rotational angular velocity are shown in Fig.[TO]as func- 
tions of the binary mass ratio, for a Roche-lobe filling 0.3Mq 
helium white dwarf in a circular binary. The different curves 
in the bottom panel of the figure correspond to white dwarf 
rotation periods ranging from 1 to lOOOhr The tidal evo- 
lution time scales for the orbital semi-major axis are orders 
of magnitude larger than those due to gravitational radiation, 
so that the effects of quasi-static tides on gravitational wave 
signals from white dwarf binaries can be neglected even for 
semi-detached white dwarf binary systems such as AM CVn 
binaries. The tidal spin-up time scale, on the other hand, can 
be shorter than 1 Myr, provided that the white dwarf is rotat- 
ing much slower than the orbital motion of the companion. 

8. CONCLUDING REMARKS 

We derived a formahsm to study dissipative tides in white 
dwarfs in the limiting case of quasi-static tides. The limit 
is applicable to binaries of arbitrary eccentricities as long as 
the forcing angular frequencies of the non-negligible terms 
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in the Fourier expansion of the tide-generating potential are 
smaller than the inverse of the white dwarf's dynamical time 
scale. We account for both convective and radiative damp- 
ing of quasi-static tides, but find the total perturbation of the 
gravitational potential, which determines the binary's tidal 
evolution, to be affected mainly by convective damping. At 
the order of approximation considered, radiative damping af- 
fects the total perturbation of the gravitational potential only 
through the surface boundary condition expressing the conti- 
nuity of the gravitational potential and its gradient at the star's 
surface. 

We applied the formalism to binaries consisting of a 0.3 Mq 
helium white dwarf and a 0.3Mq or 10^ M© point-mass com- 
panion, representative of double degenerates (e.g. Nelemans 
et al. 2001) and coalescing white dwarf-massive black hole 
binaries (Menou, Haiman, & Kocsis 2008; Sesana et al. 
2009), respectively. Since the evolutionary time scales depend 
strongly on the ratio of the white dwarf radius to the orbital 
semi-major axis, we focused the application on short-period 
binaries with orbital frequencies in LISA gravitational wave 
frequency band (10~'*-10~' Hz). The time scales for the rate 
of secular change of the orbital semi-major axis and eccen- 
tricity are found to be longer than a Hubble time for all binary 
configurations considered. Orbital evolution due to convec- 
tive damping of quasi-static tides can therefore be neglected 
in the construction of gravitational wave templates of white 
dwarf binaries for LISA. The time scales for the rate of secu- 
lar change of the white dwarf's rotational angular velocity, on 
the other hand, can be shorter than lOMyr, especially if the 
white dwarf is initially rotating with a frequency that is much 



smaller than the binary orbital frequency. 

Even though tidal spin-up of the white dwarf can occur in 
astrophysically interesting time spans, the time scales are still 
much longer than the orbital inspiral time scale due to gravi- 
tational radiation. Tides will therefore not be able to catch up 
with gravitational-radiation driven orbital evolution and syn- 
chronize the white dwarf's rotation with the orbital motion. 
White dwarf binaries in the LISA band therefore naturally 
evolve from a low- to a high-frequency tidal forcing regime. 
For the shortest period binaries considered, the quasi-static 
tide approximation therefore breaks down and energy dissi- 
pation through dynamic tides and, in particular, resonantly 
excited ^-modes must be taken into account. Possible mech- 
anisms contributing to damping of nonradial ^-modes are ra- 
diative heat leakage, neutrino losses, and gravitational radia- 
tion (Osaki & Hansen 1973; Rathore, Blandford, & Broderick 
2005), or, in the case of large amplitude modes, nonlinear cou- 
pling to other nonradial oscillation modes (e.g. Dziembowski 
1982; van Hoolst 1994ab, Wu & Goldreich 2001). We will 
explore energy dissipation through dynamic tides in detail in 
future investigations. 
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APPENDIX 

A. PERTURBATION OF THE TEMPERATURE 



Assuming the mean molecular weight of a mass element does not change on time scales comparable to or shorter than a tidal 
oscillation period, the Lagrangian perturbation of the entropy S due to a star's tidal distortion can be related to the Lagrangian 
perturbations of the temperature T and mass density p as (Unno et al. 1989, Eq. 13.74) 



{SS)nde = Cv 



iST)n 



(Sp) 



tide 



(Al) 



where Cy is the specific heat per unit mass at constant density, and a i5 in front of a quantity denotes the Lagrangian perturbation 
of that quantity. By taking the total time derivative of this expression, using the property that the total time derivative commutes 
with the Lagrangian perturbation operator, and substituting dS/dt = (1 /T){dQ/dt), it follows that 



am 


= Cv 


d 






tide 


dt 





-(r3-i) 



dt 



tide. 



(A2) 



Next, passing on from Lagrangian to Eulerian perturbations and taking into account that the unperturbed star is static and in 
thermal equiUbrium yields 



dQ 
dt 



tide 



dt \ T T dr ^"^ 



f dp 
ot \ p p dr 



(A3) 



Finally, after separating the time t and angular coordinates 9 and (j) by expanding the tidal displacement field and the perturbed 
stellar structure quantities in Fourier series similar to those used in § [3] and eliminating the Eulerian perturbation of the mass 
density by means of Eq. (fT9b , the radial part of the Eulerian perturbation of the temperature associated with the spherical harmonic 
Y^"{9, (j)) and the forcing angular frequency a,„,k is given by 



^ 



1 dT 



£,i.m.k - (J^i - 1) a£,m,A- " i 



1 



<ym.kCv T 



dQ 
dt 



(A4) 



l.m,k 



In analogy with Eq. ( |26] |, Eq. (IA4l i can be cast in dimensionless form by expressing the physical quantities in the units listed 
in Table |2] In addition, the temperature T and specific heat Cy are expressed in the units GMi /(T^gas^i) and 7?.gas, where 7?.gas is 
the universal gas constant. It follows that 



T' 

'-l.m.k 



~r~ t,i.m.k " 

T dr 



(Fa - 1)q!£,„ jt-icr,„.<: 



dQ 
dt 



(A5) 



e,in,k 
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where C is defined by Eq. ( |28] |. At 0{a^^ g^), the Eulerian perturbation of the temperature is then given by 



7^/(0) , ,^ 



and at Oicr,,,^!,) by 



T dr^'-'-* ^ ' ^ CvT\dt ) 



^— &-(r3-l)<L-iC'^ ( ] . (A7) 

l,m,k 



B. PERTURBATION OF THE RATE OF CHANGE OF THERMAL ENERGY 

The appHcation of boundary Condition ( |40] i for the total perturbation of the gravitational potential and its first derivative 
requires the evaluation of the Eulerian perturbation of the rate of change of thermal energy at the star's surface. Denoting the 
rate of energy generation per unit mass by e and the energy flux by F, the Eulerian perturbation of the rate of change of thermal 
energy is given by 

'fY =4de+%v.F-iv-f:;,„ (Bi) 

dt J tide P P 

where 

Aac 

F = Vr. (B2) 

3 Kp 

In the latter equation, a is the radiation constant, c the speed of light, T the temperature, and k the opacity. By using Eq. ( IB2I ). 
we implicitly assume the energy transfer in the near the star's surface to be radiative and/or conductive (e.g. Hansen, Kawaler, 
& Trimble 2004). Since a proper treatment of the contribution of convection to the perturbation of the rate of change of thermal 
energy requires a non-local time-dependent theory of convection, which is beyond the scope of this investigation, we here simply 
neglect convective effects in Eq. ( IB lb for {dQ/dt)[^^^. In the following derivation, all quantities are assumed to be expressed in 
the same units as those adopted in §|4](see Table|2]i. 

The Eulerian perturbation of the energy generation rate can be expressed in terms of the Eulerian perturbations of the mass 
density p, the temperature T, and the mean molecular weight p as 

£tide _ ^ Ptide , ^ -^tide , ^ P-tide tVi^\ 

e p T p 

where 

fdlne\ fdlne\ f dlne\ 

^"^'^J./ '^H^J,/ ^'^H^J,/ 

If diffuse mixing is assumed to be negligible on time scales of the order of the tidal oscillation period or shorter, the Lagrangian 
perturbation of the mean molecular weight is zero, so that 

Mt'ide=-^6ide (B5) 

dr 

(e.g. Savonije & Papaloizou 1984, Witte & savonije 1999). 

From Eq. ( IB2b and the spherical symmetry of the unperturbed star, it follows that the Eulerian perturbation of the energy flux 
is given by 



-I r 



T K p 



(B6) 



In this equation, the Eulerian perturbation of the opacity can be expressed in terms of the Eulerian perturbations of the mass 
density p, the temperature T, and the mean molecular weight p as 

^tide _ Ptide , ,^ ^tide , Ptide /■d'J\ 
K p 1 p 



where 

'dlnK\ / 51n 



After separating the time t and the angular coordinates 6 and cf) by expanding the tidal displacement field and the perturbed 
stellar structure quantities in Fourier series similar to those used in § |3] and using Eqs. (|B5b - (|B7| i, it follows that 
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Next, introducing expansions of the form given by Eq. ( |29] l, substituting the solutions for p'*"' and given by Eqs. dSTT l 
and ( |A6t , and using the chain rule yields 



dt J 



'(0) 

l.mM 

1 d 



ds 
dr 



f(0) 
'ie.m.k 



1 d 



pr^ dr 
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This expression can be simplified by noting that 
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After some algebra and using the property that dQ/dt = in the unperturbed star, it follows that 
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By using differential Eq. (|32] | and the Relation (ISTT i between '^fl, f, and ^^^^j^ j,, this expression can be further simplified to 



dQ\ 
dt J 



'(0) 



i.m.k 



F dr g dr 



F { \ dF \dg\ dS.fla 
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or, equivalently. 
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(B14) 



(B15) 



where L,- = 4tt r- F is the luminosity at radius r. Substitution of this expression into Eq. ( I38b and differentiation of the radial 
component of the tidal displacement field with respect to time yields an expression for the radial component of the tidal velocity 
field in agreement with Eq. (20) from Campbell (1984). 

C. SOLUTION OF THE SYSTEM OF DIFFERENTIAL EQUATIONS GOVERNING QUASI-STATIC TIDES 

In the quasi-static approximation, the total perturbation of the gravitational potential associated with the spherical harmonic 
Yi'{9, (f)) and forcing angular frequency cr,,, a in Expansion (O of the tide-generating potential is given by 



'£,m,A:W-riym,t e^„,^i,V J, (CI) 

where the functions ip^'^^r) and ip^^\r) are determined by differential Eqs. (|32] i and (39[ . respectively. A general solution to 

'jf\r) and V'f 

= Ci V'f V) + C2 <V), (C2) 



Eq. ( |32] | consists of a linear combination of two independent particular solutions ipf \r) and ipf \ry. 



where Ci and C2 are two undetermined constants. We let ip\ (r) be the particular solution that behaves as r^ near r = 0, and tpf'ir) 
the particular solution that behaves as r"^"' near r = 0. In order for the total perturbation of the gravitational potential and the 
radial component of the tidal displacement field to remain finite at r = 0, we set C2 = 0. The constant Ci is determined by means 
of boundary Condition (l33Tl"*: 



Cl =-etide(2^+l)Q., 



,m.k 



dip 



(0)' 



dr 



+(£+i-^)vr(i) 



(C3) 



where p.s and gs are the mass density and gravity at the star's surface. 

* Note that we use the same dimensionless quantities as those adopted in §|4] 
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Next, a general solution to Eq. ( |39l ) is given by 



A(r) 



A(r) 



(C4) 



where C[ and C2 are two undetermined constants, V'l'^'") and are two independent particular solutions of the homogeneous 

part of Eq. i39[ . and the functions (3{r) and A(r) are defined as 



1 d 
gr^ dr 



pr V 



d /CitAi 



(0)' 



dr 



and A(r) = t/i 



' t/r ^2 ■ 



(C5) 



The homogeneous part of Eq. (|39] l is formally the same as Eq 



We therefore adopt the same two independent particular 



solutions as before, and let V'j'V) be the particular solution that behaves as / near r = 0, and ?/'2''('") the particular solution that 
behaves as r"^"' near r = 0. 

From a numerical point of view, it is convenient to avoid the numerical calculation of the derivatives of the mass density p 
and turbulent viscosity coefficient v in the function /3(r) appearing in the solution for ^/'^''(r). By performing an integration by 
parts and using that t/'I' V) and V2' V) 

are solutions to the homogeneous part of Eq. ( [39] l. we therefore transform the solution for 

into the form 
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Because of the requirement that the total perturbation of the gravitational potential and the radial component of the tidal displace- 
ment field must remain finite at r = 0, we set Cj = 0. Boundary Condition ( |40] i then yields 
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where A'^^ is the square of the Brunt- Vaisala frequency at the surface of the unperturbed star, and 
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The constant C[ is thus purely imaginary. 



D. ALTERNATE FORM OF THE TIDAL EVOLUTION EQUATIONS 

Equations ( l54l i. dSSl l, and (|67] | for the rates of secular change of the orbital semi-major axis a, the orbital eccentricity e, and 
the rotational angular velocity fli are of a different form than the equations derived by, e.g., Zahn (1977, 1978), Hut (1981), and 
Ruymaekers (1992). To show that the different forms are equivalent, we restrict ourselves to the dominant £ = 2 terms and use 
Kepler's third law to rewrite Eq. ( |54l i as 



= 4^,(1.,) 
^1 



2 +00 



^^H2.„i.k \F2,m,k\ sin72,m,/tG2ji.i(e)- 



(Dl) 



m=-2 k=0 



If we furthermore assume dissipative effects to be small, the phase angles ^2,m,k are small and proportional to the forcing angular 
frequencies ai„,k, so that we can set 

sin72,,„,i w-cr,„,/t7"2- (D2) 



The constant T2 has dimensions of time and is independent of m and k because the quasi-static tide solutions ^'^'^jX^i) and 
*f!iA:(^i) are divided by C£ ,„ i in Eq. ( I6OI 1 for the phase angles 7£,„,,/t. The minus sign in Eq. (ID2l i is included to take into account 
that the tides lag behind the position of the companion when Oi < n. 
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By the use of the definition of the forcing angular frequencies ami, Eq. (IDH can then be cast in the form 
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In the limiting case where all forcing angular frequencies a,n,k tend to zero, the constants Fa./n,* all tend to the classical apsidal 
motion constant ^2 (Smeyers & Willems 2001). It can then be shown numerically or through the use of Taylor series that 
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After setting E\I{GM\ T2) = tf, where tf is a characteristic tidal energy dissipation time scale, the equation for the rate of secular 
change of the orbital semi-major axis in the limiting case of weak damping and small forcing angular frequencies takes the same 
form as the equation for the rate of secular change of the orbital semi-major axis derived by Zahn (1977, 1978), Hut (1981), and 
Ruymaekers (1992)^. 

Similarly, retaining only the dominant £ = 2 terms, Eqs. ( fSST l and (l67T i can be rewritten as 
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/f(.2)-ll(i-.2)^/Vf(.^)^' 

18 n 



dill 
dt 



= 6 



GMi _ 2^1^? f 



R 



■T2q 



/f(.2)-(i-.2)^/Vf(e')- 



where 



7 n -p2\13/2 2 +00 



4 (1-e-) 
33 7 



m=-2 k=0 



2-j5 2 +00 

E E '"^^.mX \F2.m.k\ G2^i,t(e)- 



m=-2 k=0 
2 +00 



«i=-2 k=0 
2 +OC 



/f (.2) = Ui-e'f'- E E 1^-2,...! G'Z/'^- 



(D8) 



(D9) 



(DIO) 



(DID 



(D12) 



(D13) 



n=-2 k=0 



In the limiting case where all forcing angular frequencies cr„, 1 tend to zero, it can again be shown numerically or through the use 
of Taylor series that 



/<V)^.2(l + ^.^.f.^4.^ 



f^''\e^)^k2{l + -e' + -e 



3 2 ,1 4 



/iV)^A:2( 1 + 



16 



/f (e')^fc2(^l + 3e"+^< 
^ Equation jD3t differs from Eq. (9) in Hut (1981) by a factor of 2 due to a different definition of the apsidal motion constant A'2. 
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(D14) 
(D15) 
(D16) 
(D17) 
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Hence, after setting R\ / {GMi T2) = tf the equations for the rate of secular change of the orbital eccentricity and rotational angular 
velocity in the limiting case of weak damping and small forcing angular frequencies also take the same form as the equations for 
the rate of secular change of the orbital eccentricity and rotational angular velocity derived by Zahn (1977, 1978), Hut (1981), 
and Ruymaekers (1992). 
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